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Abstract 


Tins  report  suinniarizes  researcli  conducted  luring  the  past  year  whose  aim  is  to 
develop  and  test  an  accurate  and  efficient  nieiliod  for  describing  boundary-layer 
development  in  the  transition  region.  .4  low-Reynolds-number  version  of  the 
Wilcox  k  model  has  been  postulate<l  and  tested  against  experimental  data 
for  10  fully  turbulent  boundary  layers  and  for  more  than  20  transitional  rases. 
Overall  discrepancies  between  theory  and  experiment  are  smaller  than  those  for 
simpler  correlations  designed  for  use  with  algebraic  turbulence  models.  Singular 
perturbation  analysis  of  the  compressible  vi.scous  sublayer  and  defect  layer  has 
been  initiated  with  some  interesting  preliminary  results.  Some  progress  has 
been  made  toward  eliminating  the  k  —  model's  .sensitivity  to  the  freestream 
boundary  condition  on  jj. 
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Chapter  1 

Introduction 


The  purpose  of  this  research  project  for  the  past  three  years  has  been  to  develop 
analytical  and  computational  tools  suitable  for  predicting  properties  of  hyper¬ 
sonic  turbulent  flows.  Significant  progress  has  been  made  on  several  aspects 
of  the  turbulent-flow  problem,  especially  in  representing  low-Reynolds-number 
effects  and  boundary-layer  transition.  Because  of  the  importance  of  transition 
in  hypersonic  boundary  layers,  most  of  our  research  during  the  past  year  has 
focu.sed  exclusively  the  transition  problem. 

1  he  overall  objective  of  our  approach  to  the  transition  problem  is  to  use  the 
Wilcox’  k  —  ^  turbulence  model  as  the  foundation  for  studying  and  modeling 
the  transitional  flow  region.  (?onsistent  with  the  needs  of  NASA,  the  transition 
point  is  assumed  to  be  known  a  priori.  Computations  ran  thus  be  initiated  at  the 
known  transition  location  and  continued  downstream  through  the  transitional 
flow  region  and  well  into  the  fully  turbulent  region.  To  develop  the  model,  we 
have  followed  a  sequence  of  interrelate<|  steps. 

1.  We  have  drawn  from  the  extensive  work  done  by  Wilcox'” '  to  help  for¬ 
mulate  a  low-Reynolds-number  version  of  the  h  model.  While  much 
of  this  early  research  focu.sed  on  locating  the  transition  point,  some  of 
the  analysis  (most  notably  Wilcox'’)  revealed  the  way  in  which  the  model 
equations  describe  the  growth  of  the  instability  from  the  transition  point 
to  the  fully  turbulent  regime. 

2.  We  have  simulated  several  of  the  fl()ws  that  have  been  done  with  Direct 
Numerical  Simulation  (DNS)  methods.  Tlie  results  of  Mansour,  Kim  and 
.Vloiti*  have  received  immediate  attention  for  two  different  Reynolds  num¬ 
bers.  Comparing  model  predictions  with  the  DNS  results  helped  greatly 
in  developing  the  model. 

We  have  tested  the  model  against  10  two-dimensional,  fully  turbulent 
boundary  layers,  including  both  incompressible  and  compressible  cases. 
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1.  VVe  have  tested  the  model  against  all  of  the  two-dimensional  rases  pre¬ 
senter!  by  Singer,  Dinavahi.  ami  Iyer.' 

0.  We  have  done  some  analysis  specifirally  aimeil  at  eliminating  the  <•  —  ^ 
model  s  sensitivity  to  the  frt'esiream  boiintlary  condition  i^tn  a.’. 

h.  We  have  initiated  a  singnlar  perturbation  analysis  of  tin'  compressible 
turbulent  boundary  layer  to  hel[)  riuantify  effects  of  comprr'ssibility  on 
model  predictions. 

.\ddressing  Items  1  and  2  has  produced  a  low-Reynolds-number  version  of 
the  k  —  ^  model.  Wilcox^**  describes  tht'  model  and  prt'sents  applications  to 
fully  turbulent  channel  and  pipe  flow,  and  for  a  transitional,  incompn'ssible 
flat-plate  boundary  layer.  He  also  clarifies  the  manner  in  which  the  model 
predicts  transition. 

This  report  applies  the  model  to  10  turbulent  boundary-layer  test  ca.ses  and 
to  more  than  20  transitional  flows  for  which  experimental  data  are  available. 
The  applications  show  that  additional  model  development  is  needed  to  prop¬ 
erly  account  for  effects  of  compressibility.  We  have  initiated  the  new  model 
development  efforts  by  addressing  Items  I  and  a. 

Chapter  2  sununarizes  the  low-Reynolds-number  version  of  the  k~^  model. 
Chapters  0  and  4  present  results  of  our  applications.  Boundary  condition  and 
closure  coefficient  sensitivity  is  discussed  in  Chapter  a.  We  analyze  the  com¬ 
pressible  turbulent  boundary  layer  in  ( 'hapter  fi.  Chapter  7  summarizes  results 
of  the  research. 


Chapter  2 


Low-Reynolds-Number  k-w 
Model 


The  equations  of  motion  for  the  Wilcox^®  low-Reynolds-number  k  —  ui  model 
were  initially  postulated  for  incompressible  flows.  To  generalize  for  compressible 
flows,  we  introduce  Favre^^  mass  averaged  variables  and  rewrite  the  equations 
of  motion  as  follows. 


dp  d 

-  +  gj-(„)  =  o 


d,  d 


d  p  ,  drji 


dt 


dxj 


(pUjUi)  =  -  —  + 


dxi  dxj 


(2.1) 

(2.2) 


^  {pE)  +  -^ipujH)=  ^ 


dt 


dxj 


(2.3) 


9  ,  , .  d  .  , .  r  dui  d 

jPP>‘)+—  l)  =  —  -  ,3  +  — 


dxj  L 


(/r  +  (T*/iT)7 


dxji 


(2.4) 


d  ,  .  d  ,  ,  ui  T’  dui  ,2  d 


{p  +  apr)- 


du) 


dxj\ 


(2.5) 


In  Equations  (2.1)-(2.5),  t  and  Xi  denote  time  and  position  vector;  p  and  p 
are  density  and  pressure;  Ui  is  mass-averaged  velocity  vector;  E  and  H  are  total 
energy  and  total  enthalpy;  Tij  and  (p  are  the  total  stress  tensor  and  the  total 
heat  flux  vector;  is  the  Reynolds  stress  tensor;  k  and  w  are  turbulence  kinetic 
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energy  and  specific  dissipation  rate;  and,  a,  /5,  0* ,  a,  a*  are  closure  coefficients. 
The  following  constitutive  relations  are  needed  to  close  the  system. 


♦  I  T 

PT  =  a—,  Tij=Tlj+Tf 
ul  ■' 


(2.6) 


=  2//  (  Sij  - 


1  duk 


rT.=2nT{Sij-r^^6ij)-^pk6,, 


'li  = 


_ £r  \ 

Pi’i  Prr  )  dxj 


E  =  e  +  —  UiUi  +  k,  H  =  h  +  -UiUi  +  k 

1  /  dui  duj 
H - - 


2  \dxj  dxi 


(2.7) 

(2.8) 

(2.9) 

(2.10) 

(2.11) 


In  the  constitutive  relations,  we  have  introduced  the  molecular  stress  tensor, 
rk,  mean  strain-rate  tensor,  5',j,  molecular  viscosity,  p,  eddy  viscosity,  pTi 
enthalpy,  h,  laminar  Prandtl  number,  Pri,  turbulent  Prandtl  number,  Prr, 
internal  energy,  e,  and  an  additional  closure  coefficient,  a*.  Finally,  to  complete 
closure  of  the  system,  we  specify  the  values  of  the  closure  coefficients  as  follows. 


,  «*  +  Ef‘-Tl  Ih 

1  +  Rer/Rk 

5  ao  +  Rer/Rw  ,  .,_i 
"=9'  I  +  Kct/R..  ' 

(2.12) 

9  5/lS  +  {ReT/Rp)^ 

^  100  l+iRer/Rpy 

(2.13) 

,  3 

^  40’ 

.  1  1 
"  =2’  "^2 

(2.14) 

/^/3,  a,  -  1/10 

(2.15) 

R0  =  6, 

Rk  =8,  /?u>  =  27/10 

(2.16) 

The  quantity  Rer  is  the  turbulence  Reynolds  number  defined  by 

II 

(2.17) 

•jjji 
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Chapter  3 

Turbulent  Boundary-Layer 
Applications 


The  primary  objective  of  this  research  project  is  to  describe  boundary-layer  de¬ 
velopment  through  transition  from  laminar  to  turbulent  flow.  As  part  of  this 
objective,  it  is  important  that  we  provide  an  accurate  description  in  the  tur¬ 
bulent  region  immediately  following  transition,  i.e.,  we  insist  that  our  model 
approach  the  proper  limiting  state  of  the  boundary  layer.  Consequently,  since 
the  k  —  ui  model  without  viscous  corrections  generally  is  very  accurate  for  turbu¬ 
lent  boundary  layers,  a  round  of  tests  is  in  order  to  make  sure  model  predictions 
are  not  adversely  affected  by  the  viscous  corrections.  We  have  performed  ten 
boundary  layer  computations  including  effects  of  adverse  and  favorable  pressure 
gradient,  and  for  Mach  numbers  as  high  as  10.  Table  3.1  lists  the  ten  cases. 


Table  3.1:  Turbulent  Boundary  Layer  Test  Cases 


Description 

Vp 

Data  Source 

Incompressible  Flat  Plate 

None 

Wieghardt-Tillman^^ 

Flow  1100 

Adverse 

Ludweig-Tillman^^ 

Flow  1200 

Adverse 

Ludweig-Tillman^^ 

Flow  1300 

Favorable 

Ludweig-Tillman^^ 

Bradshaw  Flow  C 

Adverse 

Bradshaw^'* 

Samuel- Joubert 

Adverse 

Samuel-Joubert^® 

Mach  2.244  Flat  Plate 

None 

Shutts^® 

Mach  4.544  Flat  Plate 

None 

Coles^® 

Mach  10.31  Flat  Plate 

None 

Watson 

Mach  2.653  Boundary  Layer 

Adverse 

Fernando-Smits^^ 
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Incompressible  Flat  Plate.  Figure  3.1  compares  computed  and  measured 
skin  friction,  Cf,  and  sublayer-scaled  velocity,  !«+.  As  shown,  the  only  significant 
difference  between  predictions- with  and  without  viscous  corrections  is  in  the 
transition  point.  The  unmodified  model  undergoes  transition  much  earlier  than 
the  model  with  viscous  corrections.  Skin  friction  at  the  final  station  is  2.50- 10“^ 
with  viscous  modifications  compared  to  2.49  •  10“'^  without. 

Flows  1100  and  1200.  Figures  3.2  and  3.3  compare  computed  and  mea- 
sured*^^  Cf  and  for  two  incompressible  boundary  layers  with  adverse  pressure 
gradient,  flow  1100  has  a  mild  adverse  gradient  while  Flow  1200  has  a  fairly 
strong  adverse  gradient.  Note  that  computations  have  been  terminated  at  the 
last  data  station  upstream  of  a;  =  3  m.  for  both  flows.  As  noted  by  Coles  and 
Hirst, the  data  fail  to  satisfy  the  momentum  integral  equation  downstream  of 
X  =  3  m.  and  are  thus  unreliable.  As  with  the  flat  plate  case,  differences  between 
results  obtained  with  and  without  viscous  corrections  are  barely  noticeable.  For 
Flow  1100,  skin  friction  at  the  final  station  is  1 .90  •  10“''^  with  viscous  corrections 
and  1.93  •  10~^  without.  A  slightly  larger  difference  is  found  for  Flow  1200, 
namely,  Cf  —  1.58  •  10“'^  with  viscous  corrections  and  1.61  ■  10“^  without. 

Flow  1300.  Figure  3.4  compares  computed  and  measured  properties  for 
Flow  1300,  another  one  of  the  incompre,ssible  Ludweig-Tillman^^  flows,  this 
time  with  favorable  pressure  gradient.  Velocity  profile  differences  are  noticeable 
only  below  y"*"  =  200,  and  skin  friction  is  nearly  identical.  Skin  friction  at  the 
final  station  is  3.34- 10"^  with  viscous  corrections  included  and  3.35- 10“^  when 
they  are  excluded. 

Bradshaw  Flow  C.  This  incompressible  adverse  pressure  gradient  case  was 
one  of  the  most  difficult-to-predict  case  in  Stanford  Olympics  As  shown 
in  Figure  3.5,  with  or  without  viscous  corrections,  the  k  —  u)  model  virtually 
duplicates  measured  cj  and  for  this  flow.  At  the  final  station,  cj  =  1.60-10“^ 
with  viscous  corrections  and  1.61  ■  10“'^  without. 

Samuel- Joubert.  This  incompressible  adverse  pressure  gradient  boundary 
layer  was  supposed  to  be  a  simple  application  in  Stanford  Olympics  11.^®  On 
the  contrary,  it  proved  to  be  very  difficult  for  all  models  and  has  become  a  key 
test  case  for  how  well  a  turbulence  model  predicts  effects  of  adverse  pressure 
gradient.  Figure  3.6  shows  that  both  high-  and  low-Reynolds-number  versions 
of  the  model  nearly  duplicate  measured  skin  friction  and  velocity  profile  at  the 
la§t  station.  At  the  final  station,  cj  —  1.21  ■  10“®  with  viscous  corrections  and 
1.25  •  10“^  without. 

Compressible  Flat  Plates.  As  shown  in  Figures  3.7,  3.8  and  3.9,  the 
viscous  corrections  have  virtually  no  effect  on  skin  friction  and  velocity  profiles 
in  the  fully  turbulent  region  for  these  three  applications.  For  all  three  cases, 
skin  friction  differs  by  less  than  one  percent  at  the  final  station.  This  contradicts 
the  claims  of  Zhang,  et  aP°  that  achieving  correct  asymptotic  near-wall  behavior 
is  required  to  accurately  predict  compressible  flows.  If  their  claims  were  correct, 
the  Ar  — w  model:  (a)  would  not  predict  flow  properties  so  close  to  measurements 
for  the  Mach  2.244  and  4.544  cases;  and,  (b)  would  be  closer  to  measurements 
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Kigure  3.1:  Incompressible  fiat-plate  boundary  layer;  o  Wieghardt-Tillman 
h  —  ^  with  viscous  corrections: - k  ~  u)  without  viscous  corrections. 


Figure  3.4:  Flow  1300:  Incompressible  boundary  layer  with  favorable  pressure 
gradient;  o  Ludweig-Tillman;  —  k  —  uj  with  viscous  corrections;  —  —  —  k  —  u> 
without  viscous  corrections. 


Figure  3.5;  Bradshaw  Flow  C:  Incompressible  boundary  layer  in  adverse  pres¬ 
sure  gradient;  o  Bradshaw;  —  k  —  u>  with  viscous  corrections;  —  —  —  k  —  u 
without  viscous  corrections. 


Figure  3.6:  Incompressible  boundary  layer  in  an  increasingly  adverse  pressure 
gradient;  o  Sainuel-Joubert;  —  k  —  ui  with  viscous  corrections;  —  —  —  k  —  uj 
without  viscous  corrections. 


for  the  Mach  10.31  case  witli  viscous  corrections  included.  What  their  analysis 
actually  shows  that  they  have  simply  canceled  the  k  —  t  model’s  inherent  flaw 
in  defect-layer  structure^  '^  by  projecting  their  viscous  modifications  far  above 
the  viscous  sublayer. 

Mach  2.563  Boundary  Layer.  The  final  turbulent  flow  application  is  for 
a  compressible  boundary  layer  with  a  freestream  Mach  number  of  2.563.  The 
boundary  layer  is  subjected  to  an  adverse  pressure  gradient.  Figure  3.10  shows 
that,  as  in  virtually  all  of  our  applications,  only  slight  differences  are  present 
between  model  predictions  with  and  without  viscous  corrections.  At  the  final 
station,  skin  friction  with  viscous  corrections  is  1.1110"'^  and  1.1010“^  without. 

In  summai-y,  for  all  ten  cases  considered,  differences  between  computed 
flow  properties  with  and  without  viscous  corrections  are  less  than  3%,  and  are 
generally  less  than  1%.  Hence,  the  viscous  corrections  leave  the  best  features  of 
the  k  —  oj  model  intact,  i.e.,  the  model  still  accurately  predicts  effects  of  pressure 
gradient  and  compressibility  up  to  Mach  5. 
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Figure  :?.7:  Mach  ‘2,244  adiabatic-wall  Hat-plate  boundary  layer;  o  Shutts; 
--  k  —uj  with  viscous  corrections; - k  -  u)  without  viscous  corrections. 
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Chapter  4 

Transitional 

Boundary-Layer 

Applications 


In  order  to  test  the  suitability  of  the  low-Reynolds-number  k—<j  model  for  tran¬ 
sitional  boundary-layers,  we  have  computed  all  of  the  two-dimensional  transi¬ 
tional  boundary  layer  cases  considered  by  Singer,  et  al.®  Table  4.1  summarizes 
the  cases  we  have  done. 

Table  4.1:  Transitional  Boundary  Layer  Test  Cases 


Flow 

Description 

Vp 

Data  Source 

1 

Incompressible  Flat  Plate 

None 

Schubauer-KlebanofP^ 

2 

Favorable  Vp 

Favorable 

Blair-Werle^^ 

3 

Supersonic  Cone  Flow 

None 

Fisher- Dougherty^® 

3 

Supersonic  Cone  Flow 

None 

Chen,  et  al^^ 

4 

Freestream  Turbulence 

None 

Schubauer-Skramstad^® 

4 

Freestreanj:  Turbulence 

None 

Blair^^ 

'  5 

Prolate  Spheroid 

Favorable 

Meier-Kreplin-Ming^® 

6 

Concave  Surface 

None 

Swearingen-Blackwelder^® 

As  discussed  by  Wilcox^®,  transition  predictions  with  the  k  —  u:  model  are 
sensitive  to  the  freestream  values  of  k  and  w,  most  notably  to  the  former.  All 
of  the  computations  have  been  done  using  the  following  boundary  condition  at 

•X 
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the  start  of  the  computation: 


where  S  denotes  boundary  layer  thickness,  and  subscript  e  denotes  boundary- 
layer  edge.  Because  our  goal  is  to  examine  model  predictions  in  the  transition 
region  as  opposed  to  predicting  the  actual  transition  point,  the  freestream  tur¬ 
bulence  level  has  been  selected  to  match  the  measured  transition  location  for 
most  of  the  computations.  We  quote  the  initicd  value  of  the  freestream  intensity, 
T',  for  all  cases,  where 

T'  =  lOOyf^  (4.2) 

The  quantity  (7e  is  boundary-layer  edge  velocity.  The  freestream  values  of  k 
and  ui  vary  throughout  the  computation,  and  follow  from  solution  of  the  model 
equations  in  the  freestream.  See  Wilcox'’  for  additional  details. 

In  all  cases,  computation  begins  at  the  plate  leading  edge,  and  the  turbulence 
kinetic  energy  is  initially  set  to  the  freestream  value  throughout  the  boundary 
layer.  The  initial  w  profile  is  given  by  the  exact  laminar-flow  solution  to  the 
model  equations. 

Incompressible  Flat  Plate.  The  first  case  is  Flow  1  from  Singer,  et  al.^ 
This  is  an  incompressible  flat-plate  boundary  layer  that  undergoes  transition 
at  a  plate  length  between  1.6  m.  and  1.8  m.  The  initial  freestream  intensity  is 
assumed  to  be  0.043%,  and  decreases  with  distance  along  the  plate  since  the 
solution  to  the  model  equations  in  the  freestream  tells  us  the  turbulence  decays. 
At  the  end  of  the  computation  .the  freestream  intensity  has  dropped  to  0.025%. 
The  nominal  value  quoted  for  the  experiment  is  0.030%. 

Figure  4.1  compares  computed  and  measured  skin  friction  throughout  the 
transition  region.  Computed  and  measured  c/  differ  by  less  than  14%  of  the  peak 
skin  friction.  According  to  Singer,  et  al,  the  Dey-Narasimha''®  correlation  yields 
values  about  6%  lower  than  measured  for  this  flow  while  the  ONERA/CERT^® 
model  predicts  peak  skin  friction  about  25%  higher  than  measured. 

Favorable  Pressure  Gradient  Boundary  Layers.  The  next  applications 
are  for  incompressible  boundary  layers  in  a  favorable  pressure  gradient.  The 
bofundary  layers  considered  correspond  to  Flow  2/Cases  1  and  2  of  Singer,  et  al. 
In  addition  to  having  adverse  pressure  gradient,  the  surface  is  cooled.  Figure  4.2 
compares  computed  and  measured  Stanton  number  for  the  two  cases.  Case  2 
has  a  stronger  favorable  pressure  gradient  than  Case  1.  The  initial  freestream 
turbulence  intensities  required  to  match  the  measured  transition  point  are  0.61% 
and  0.91%,  respectively,  for  Cases  I  and  2. 

Singer,  et  al  show  that  the  Dey-Narasimha  and  the  ONERA/CERT  models 
are  almost  as  close  to  the  data  as  the  k  prediction,  although  both  predict 
a  more  rapid  approach  to  the  turbulent  state.  By , contrast,  for  Case  2  the  Dey- 
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Figure  4,1:  Transitional  incompressible  flat-plate  boundary  layer. 


Narasirnha  model  predicts  a  peak  Stanton  number  nearly  double  the  measured 
value.  Also,  the  ONERA/CERT  model  fails  to  predict  transition  for  Case  2. 

Both  of  these  flows  have  a  very  high  measured  freestream  turbulence  level 
(2.0%),  which  probably  accounts  for  the  departure  of  the  skin  friction  from  the 
laminar  value  ahead  of  transition.  This  is  probably  buffeting  and  it  tends  to 
smear  the  Stanton  number  distribution.  While  the  k  —  u)  model  prediction  fails 
to  duplicate  this  feature,  the  computation  predicts  an  asymptote  to  the  fully 
turbulent  value  in  a  sensible  distance.  By  contrast,  the  Dey-Narasimha  model 
does  not  approach  the  data  until  about  1.5  m. 

Supersonic  Cone  Flow.  For  this  application  we  consider  Flow  3  of  the 
Singer,  et  al  study.  Cases  1-4  focus  on  flow  past  a  5°  half  angle  cone  with  Mach 
nutnbers  ranging  from  1.16  to  1.86,  corresponding  to  measurements  of  Fisher 
and  Dougherty.^®  Mpasurements  are  available  only  for  the  beginning  and  end 
of  transition.  Table "4.2  summarizes  the  freestream  turbulence  intensity  used, 
transition  Reynolds  number,  Rcx,,  and  Reynolds  number  based  on  transition 
width,  7?e^a,,. 

Cases  5-7  are  for  Mach  3.36  flow  past  a  5“  half-angle  cone,  with  an  adia¬ 
batic  surface.  The  three  cases  considered  have  different  unit  Reynolds  numbers. 
Figure  4.3  compares  computed  and  measured  recovery  factor,  r,  for  the  three 
cases;  freestream  turbulence  intensity  in  the  computations  is  0.1%  for  all  three 
cases.  Note  that  unit  Reynolds  number  is  smallest  for  Case  5  and  largest  for 
Case  7.  Interestingly,  consistent  with  measurements,  the  asymptotic  value  of 
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Table  4.2:  Results  for  Fisher-Dougherty  Test  Cases 


Case 

T 

Rext 

)eirp 

1 

1.16 

.0020% 

7.02  •  10“ 

0.89  ■  10“ 

0.69  •  10“ 

2 

1.30 

.0075% 

5.61  ■  10“ 

0.84  ■  10“ 

0.97  •  10® 

3 

1.55 

,0027% 

7,89  •  10“ 

1.08  ■  10® 

1.24  •  10“ 

4 

1.86 

.0073% 

7.29  ■  10“ 

1.07  ■  10“ 

1..56  ■  10“ 

the  recovery  factor  is  about  0,875  for  all  three  cases  compared  to  the  generally 
accepted  value  of  0,89.  By  contrast,  the  Dey-Narasimha  and  ONERA/CERT 
models  predict  an  asymptotic  value  in  excess  of  0.88.  As  a  simple  numerical 
experiment,  we  have  run  these  computations  without  transverse  curvature,  i.e., 
as  two-dimensional  flows.  The  asymptotic  value  of  the  recovery  factor  is  0,89. 

Freestream  Turbulence  Effects.  We  turn  now  to  effects  of  freestream 
turbulence,  i.e..  Flow  4  of  the  Singer,  et  al  stndy.  All  of  the  cases  considered 
are  for  incompressible  boundary  layers.  Table  4.3  shows  that  computed  tran¬ 
sition  width  for  Cases  1-5  is  consistently  smaller  than  measured,  especially  at 
the  smallest  transition  Reynolds  numbers.  Note  that  the  measured  freestream 
turbulence  intensity  has  been  used  in  the  compntations. 


Table  4.3:  Results  for  Schubauer-Skramstad  Test  Cases 


Case 

T' 

Rext 

Xt 

Aj:*  )ea;p 

1 

.042% 

2.85  • 

10“ 

0.88  • 

10“ 

1.00- 

■To^ 

2 

,100% 

2.75- 

10® 

0.86- 

10“ 

1.20- 

10“ 

3 

.200% 

2.20  ■ 

10“ 

0.78  ■ 

10“ 

1.50- 

10“ 

4 

.260% 

1.80  ■ 

10® 

0.69  • 

10“ 

1.40  • 

10“ 

5 

.340% 

1.40- 

10“ 

0.58  ■ 

10“ 

1.25  ■ 

10“ 

For  Cases  6,  7  and  8,  Figure  4.4  shows  that  predicted  peak  skin  friction 
is  pbout  10%-15%  lower  than  measured.  The  ONERA/CERT  model’s  peak 
skin  friction  is  about  5%-10%  higher  than  measured,  while  the  Dey-Narasimha 
model  yields  peak  skin  friction  that  is  20%-25%  higher  than  measured.  As 
with  Cases  1-5,  the  computations  match  the  measured  freestream  turbulence 
intensity,  which  is  0.25%,  1.00%  and  2.00%  for  Cases  6,  7  and  8,  respectively. 

Incompressible  Flow  Past  a  Spheroid.  This  case  is  Flow  5  from  Singer, 
et  al.  The  flow  examines  the  effect  of  a  roughness  strip  on  the  transition  of  a 
boundary  layer  over  a  prolate  spheroid  at  zero  angle  of  attack.  Experimental 
data  have  been  provided  by  Meier,  Kreplin  and  Table  4.4  summarizes 

the  four  cases,  including  the  value  of  T'  used  in  the  computations. 
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Figure  4.3:  TYansitional  flow  on  a  5°  half-angle,  adiabatic  cone  at  Mach  3.36. 
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Figure  4.4:  Transition  width  for  fiat-plate  boundary  layers. 
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Figures  4.5  through  4.8  compare  computed  and  measured  skin  friction  (based 
on  local  boundary-layer  edge  velocity)  as  a  function  of  arc  length  along  the 
spheroid.  As  shown,  computed  and  measured  skin  friction  are  closest  when 
transition  is  induced  by  the  roughness  strip.  Without  the  roughness  strip,  the 
predicted  transition  occurs  more  abruptly  than  measured. 


Table  4.4:  Meier- Kreplin- Ming  Test  Cases 


Case 

Description 

T 

1 

Uoo  =  20  m/sec, 

no  roughness  strip 

.017% 

2 

Uoo  =  20  m/sec, 

with  roughness  strip 

.153% 

3 

Uoo  =  30  m/sec, 

no  roughness  strip 

.010% 

4 

Uoo  =  30  m/sec. 

with  roughness  strip 

.115% 

Concave  Surface  Boundary  Layer  The  final  application  is  for  incom¬ 
pressible  flow  over  a  concave  surface.  Because  the  physical  flow  has  longitu¬ 
dinal  vortices,  this  application  would  be  more  appropriately  done  in  a  three- 
dimensional  computation.  Nevertheless,  following  Singer,  et  al,  we  present  re¬ 
sults  of  a  two-dimensional  computation  in  Figure  4.9.  Freestream  turbulence 
intensity  in  the  computation  is  1.21%.  The  strongest  statement  we  can  make 
is  that  the  peak  value  of  cj  matches  the  measured  value.  The  Dey-Narasimha 
and  ONERA/CERT  models  yield  peak  values  about  6%  and  16%  higher  than 
measured,  respectively. 

In  summary,  while  the  differences  between  computed  and  measured  skin 
friction  and  Stanton  number  are  less  than  those  with  the  Dey-Narasimha  and 
ONERA/CERT  models,  the  formulation  is  imperfect.  This  is  most  evident  for 
the  cases  with  large  transition  Reynolds  numbers.  In  these  cases,  the  model- 
predicted  transition  occurs  more  abruptly  than  measured,  especially  for  the 
compressible  cases.  Chapter  5  explains  why  this  is  true  and  shows  that  a  better 
choice  of  the  low-Reynolds-number  closure  coefficients  Rj,  and  may  reconcile 
the  differences. 
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Figure  4.5:  Flow  past  a  prolate  spheroid;  =  20  m/sec;  no  roughness  strip. 
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Figure  4.6;  Flow  past  a  prolate  spheroid;  U^a  —  20  m/sec;  with  roughness  strip. 
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Figure  4.7:  Flow  past  a  prolate  spheroid;  Uco  —  30  m/sec;  no  roughness  strip. 


ID^Cf 

D 


Figure  4,8;  Flow  past  a  prolate  spheroid;  Uoo  ~  30  m/sec;  with  roughness  strip. 


Figure  4.9:  'irrausitional  How  on  a  concave  surface;  incompressible  flow. 
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Chapter  5 

Boundary  Condition  and 
Closure  Coefficient 
Sensitivity 


In  this  chapter,  we  investigate  the  sensitivity  of  fc  —  w  model  solutions  to  the 
freestream  value  of  w.  The  sensitivity  is  related  to  the  physical  time  scales  rel¬ 
evant  in  shear  flows.  A  solution  similar  to  Menter’s  postulate  has  been  tested. 
Additionally,  we  have  determined  the  sensitivity  of  model- predicted  transition 
width  to  the  low-Reynolds-number  closure  coelEcients.  We  have  found  that 
transition  width  is  insensitive  to  three  of  the  postulated  low-Reynolds-number 
closure  coefficients  («„,  a*  and-  Rp)  and  is  strongly  affected  by  two  others  (Rj 
and  Ru,).  As  a  result  of  the  understanding  gained,  we  have  found  an  ad  hoc  mod¬ 
ification  that  yields  more  realistic  transition  width  for  high  transition  Reynolds 
numbers. 


5.1  Boundary  Condition  Sensitivity 

Unlike  other  two-eqjiation  models,  the  k  —  u  model  reflects  a  sensitivity  to 
freestream  turbulence.  This  sensitivity  is  physically  realistic  for  boundary  layers 
as  turbulent  boundary- layer  skin  friction  is  known  to  increase  with  increasing 
freestream  turbulence  level.  The  model  does  indeed  predict  an  increase  in  skin 
friction.  However,  the  model  also  predicts  a  strong  effect  of  the  freestream  value 
of  w  on  the  spreading  rate  of  mixing  layers,  jets  and  wakes.  Measurements  do 
not  suggest  such  sensitivity.  FinzJly,  the  sensitivity  mitigates  one  of  the  k  —  ui 
model’s  strongest  assets  in  transition  prediction,  viz,  the  ease  with  which  initial 
profiles  can  be  established.  While  the  w  equation  has  a  well-defined  laminar 
flow  solution,  the  solution  is  somewhat  sensitive  taThe  freestream  value  of  w. 
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Hence,  transition  predictions  are  affected  by  the  model’s  sensitivity  to  freestream 
boundary  conditions. 

Recall  that  for  fully  turbulent  flows,  the  high-Reynolds-number  A  —  w  model 
postulates  that  the  Reynolds  shear  stress  is  given  by 


—  v/v'  = 


k  du 
ijj  dy 


(5.1) 


As  we  approach  the  freestream,  we  expect  on  physical  grounds  that  k  and  u'v' 
approach  negligibly  small  values  for  a  fully  turbulent  boundary  layer.  It  is 
possible  to  achieve  this  according  to  Equation  (5.1)  even  with  a  small  freestream 
value  of  10,  provided  the  ratio  of  1:  to  w  remains  small. 

Having  the  ratio  of  1:  to  w  small  is  not  the  only  physical  consideration  regard¬ 
ing  a  possible  lower  bound  on  the  freestream  value  of  w.  To  understand  why, 
note  that  the  quantity  w  is  the  reciprocal  of  the  time  it  takes  the  turbulence  to 
adjust  to  changes  in  its  environment.  By  contrast,  the  mean  velocity  gradient, 
du/dy,  is  the  reciprocal  of  the  time  scale  on  which  the  mean  flow  is  changing. 
The  Reynolds  averaging  procedure  used  to  arrive  at  our  model  equations  im¬ 
plicitly  assumes  that  mean-flow  properties  are  not  changing  more  rapidly  than 
the  turbulent  fluctuations.  Hence,  we  have  made  an  implicit  assumption  that 
in  a  shear  flow. 


^min  ^ 


du 

dy 


(5.2) 


Menter^^  has  shown  that  by  computing  the  Reynolds  shear  stress  according 


to 


u'v'  =  mm 


k  O.ik  du 
Lo  ’  du/dy  dy 


(5.3) 


k—w  model  predictions  can  be  brought  into  closer  agreement  with  measurements 
for  some  flows.  This  is  tantamount  to  saying  that 


_  10  du 
^min  “  *0 

3  dy 

when  computing  thd'  eddy  viscosity. 

'  Of  course,  Menter’s  correction  was  designed  to  improve  the  model’s  predic¬ 
tions  for  separated  flows,  rather  than  to  address  boundary  condition  sensitivity. 
Nevertheless,  we  have  considered  the  possibility  that  Menter’s  correction  may 
also  remove  some  of  the  sensitivity.  To  do  so,  we  have  tested  Menter’s  proposal 
for  flat-plate  boundary  layer  flow  and  the  turbulent  far  wake.  We  have  found 
that  the  Menter  correction  does  very  little  to  remove  the  model’s  sensitivity  to 
the  freestream  value  of  w  for  either  flow. 

While  it  is  easy  to  implement,  this  approach  is  not  entirely  self  consistent. 
We  are  still  permitting  u  to  assume  values  less  tharf  ^min  and  simply  adjusting 
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the  eddy  viscosity  according  to  Equation  (5.3).  To  eliminate  the  inconsistency, 
we  can  define  a  modified  specific  dissipation  rate,  w,  according  to: 

(5.5) 


0/  =  max 


,  du 
w,  A— 
dy 


where  A  is  a  coefficient  to  be  determined.  The  eddy  viscosity,  k  equation 
equation  then  become 


(It  =  a* 


pk 

U) 


and  uj 
(5.6) 


^  ^  /  -1  T 


T 


dxj 


■  PpQu  + 


A 

dxi 


/  * 


dcj 

(p  +  apT)^^ 


(5.7) 

(5.8) 


Using  A  <  2,  this  modification  has  virtually  no  effect  on  solutions.  However, 
when  A  >  3,  boundary-layer  solutions  exhibit  no  sensitivity  to  the  freestream 
value  of  w.  However,  this  prescription  is  completely  ineffective  for  the  far 
wake.  Consequently,  we  are  still  seeking  a  satisfactory  resolution  of  the  model’s 
freestream  boundary  condition  sensitivity. 


5.2  Closure  Coefficient  Sensitivity 


Turning  now  to  closure-coefficient  sensitivity,  we  have  found  that  when  transi¬ 
tion  occurs  at  relatively  low  Reynolds  numbers  (i.e.,  typically  for  <  10®) 
the  predicted  extent  of  the  transition  region  is  consistent  with  measurements. 
By  contrast,  at  higher  transition  Reynolds  numbers,  predicted  transition  width 
generally  is  smaller  than  measured.  This  problem  is  especially  noticeable  for 
compressible  boundary  layers.  For  example,  Figure  5.1  compares  computed  and 
measured  Stanton  number,  St,  for  Mach  20  flow^^  past  a  5°  half-angle  cone. 
As  shown,  the  computed  Stanton  number  increases  much  more  abruptly  than 
measured. 

To  understand  why  the  low-ReynoIds- number  version  of  the  k—u  model  un¬ 
derpredicts  transition  width  at  high  Reynolds  number,  recall  that  the  boundary- 
la^er  form  of  the  equations  for  k  and  ui  are: 


dk  dk  f du\^  1  ^  \/  * 

du  du  w  ^  ^  2  ^  \ /  ^ 


pT  —  Oi 


(5.9) 

(5.10) 

(5.11) 
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Figure  5.1;  Mach  20  flow  past  a  5°  half-angle  cone. 


Additionally,  the  closure  coefficients  a* ,  a,  j3* ,  /?,  rr*  and  <t  are  given  by  Equa¬ 
tions  (2.12-2.17).  Inspection  of  Equations  (5.9)-(5.11)  shows  that  the  five  low 
Ileynolds  number  specific  closure  coefficients,  ocg,  «*,  Rfj,  ttk  ^nd  /i^l  accomplish 
two  ends.  First,  a*  and  «o  control  the  points  at  which  k  and  w  are  amplified. 
The  coefficients  Rk  and  control  the  rate  at  which  a,  a*  and  approach 
their  fully  turbulent  values. 

Closer  examination  of  solutions  to  Equations  (5.9)-(5.11)  reveals  that  the 
primary  role  of  R^  is  to  establish  the  peak  value  of  k  in  the  fully  turbulent 
region.  However,  it  plays  a  more  or  less  passive  role  in  the  transition  region. 
Decreasing  the  value  of  tends  to  increase  the  predicted  width  of  the  transition 
zone  slightly. 

The  value  of  a*  has  been  chosen  to  guarantee  that  small  disturbances  begin 
to  be  amplified  at  the  minimum  critical  Reynolds  number  for  the  Blasius  bound¬ 
ary  layer.  At  low  Reynolds  numbers,  the  value  of  determines  the  width  of  the 
trinsition  zone.  To  understand  this,  recall  from  Wilcox^^  that  for  the  Blasius 
boundary  layer,  the  production  terms  in  the  k  and  w  equations  first  exceed  their 
corresponding  dissipation  terms  at  the  following  “critical”  Reynolds  numbers. 


(fie®)*  — 


8100 


12150 

'  /■ 


(5.12) 

(5.13) 


33 


If  we  approximate  that  the  dominant  meciianism  in  the  early  stages  of  tran¬ 
sition  is  the  balance  amongst  the  streamwise  convection,  production  aind  dissi¬ 
pation  terms,  the  k  and  w  equations  simplify  to 


^  dx  3 


3<  l^u/dy^l  ^ 


kij) 


du 

«—  =  0 

ax 


5  / du/dy\' 

^  0  \  w  ) 


ui 


(5.14) 

(5.15) 


Since  the  ratio  of  du/dy  to  w  increases  as  for  the  Blasins  boundary 

layer,  the  following  observations  can  be  made.  First,  provided  a*  and  a,,  are 
chosen  such  that  k  is  amplified  before  u>  is  amplified,  k  will  grow  with  increasing 
values  of  x  beyond  {Rex)k-  Second,  the  V£due  of  ag  controls  the  location  at 
which  UI  begins  to  grow,  which  ultimately  terminates  the  growth  of  k,  i.e., 
determines  the  end  of  transition.  Third,  when  the  transition  Reynolds  number 
is  of  the  order  of  1  ■  10®,  the  difference  between  {Rex)k  and  {R£x)ui  is  roughly 
7%  of  the  measured  width  of  the  transition  region.  Fourth,  as  the  transition 
Reynolds  number  increases,  the  difference  between  (Re^Jk  and  (R£i;)w  becomes 
a  much  smaller  fraction  of  the  transition  width. 

Hence,  we  should  expect  the  relative  values  of  a*  and  to  have  little 
effect  on  the  width  of  the  transition  zone  at  high  Reynolds  numbers.  Numerical 
experimentation  confirms  that  at  high  Reynolds  numbers,  transition  width  is 
essentially  independent  of  the  values  of  a*  and  Og. 

We  have  also  examined  the  sensitivity  of  transition  width  to  the  two  re¬ 
maining  low  Reynolds  number  closure  coefficients,  Rk  and  R^.  In  contrast  to 
the  relative  insensitivity  to  values  of  the  other  three  low  Reynolds  number  pa¬ 
rameters,  transition  width  is  very  strongly  affected  by  Rk  and  Rgj.  Figure  5.2 
illustrates  the  effect  of  tripling  the  values  of  Rk  and  for  the  Mach  20  flow 
discussed  above.  Interestingly,  not  only  is  the  transition  width  more  accurately 
represented,  so  is  the  heat  transfer  in  the  fully  turbulent  region. 

This  exercise,  although  ad  hoc  in  nature,  points  to  a  possible  resolution  of  the 
model’s  inaccuracy  for  transition  width  at  high  Reynolds  numbers.  There  are 
indeed  pairs  of  valued  of  Rk  and  R^,  other  than  those  quoted  in  Equation  (2.16), 
that  yield  satisfactory  results  for  general  boundary-layer  applications.  Alterna¬ 
tive  values  follow  from  a  perturbation  analysis  of  the  viscous  sublayer.  Such 
analysis  is  the  topic  of  the  next  chapter. 
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Chapter  6 

Singular  Perturbation 
Analysis  of  the 
Compressible  Turbulent 
Boundary  Layer 


As  shown  in  Chapter  4,  the  low-Reynolds-number  k  -  u  model  yields  transi¬ 
tion  predictions  consistent  with  incompressible  flat-plate  boundary  layer  mea¬ 
surements,  including  sensitivity  to  freestream  turbulence  intensity.  While  the 
model  provides  a  realistic  simulation  of  the  transitional  region  at  high  freestream 
turbulence  levels  (e.g.,  >  0.5%  freestream  intensity),  the  model  predicts  more 
abrupt  than  measured  transition  for  smaller  turbulence  levels.  As  shown  m 
Chapter  5,  the  model  displays  a  strong  (possibly  nonphysical)  sensitivity  to 
freestream  turbulence  properties  for  free  shear  flows.  The  sensitivity  is  present 
for  boundary  layers  as  well,  although  the  effect  is  not  as  pronounced  (no  more 
than  10%  changes  in  skin  friction,  for  example)  as  it  is  for  free  shear  flows. 

For  all  three  cases,  viz,  transitional  bomidary  layers,  free  shear  flows  and 
turbulent  boundary  layers,  sensitivity  to  freestream  turbulence  is  felt  through 
the  entrainment  process  at  the  turbulent /nonturbulent  interface.  Specific  to  the 
transitional  case,  we  usually  speak  of  the  boundary  layer’s  recepUviiy.  Clearly, 
developing  a  satisfactory  engineering  model  for  describing  the  transition  pro¬ 
cess  requires  a  satisfactory  description  of  the  boundary  layer’s  receptivity  to 
freestream  disturbances. 

In  order  to  analyze  this  phenomenon,  the  ideal  tool  is  a  perturbation  so¬ 
lution  for  the  turbulent  boundary  layer.  Wilcox^  has  done  such  a  solution  for 
incompressible  flow,  including  effects  of  pressure  gradient.  Since  our  goal  in  this 
Contract  is  to  develop  a  model  for  compressible  flow^,-  a  compressible  defect-layer 


36 


solution  is  needed.  Hence,  we  have  initiated  development  of  such  a  solution. 

As  shown  by  Wilcox, the  perturbation  analysis  develops  an  inner  expansion 
corresponding  to  the  viscous  sublayer  and  an  outer  expansion  corresponding  to 
the  defect  layer.  In  this  chapter,  we  examine  the  viscous  sublayer  solution  first. 


6.1  The  Viscous  Sublayer 


Our  primary  motivation  for  analyzing  the  compressible  sublayer  is  to  help  quan¬ 
tify  solution  sensitivity  to  Rk  and  when  compressibility  effects  are  signifi¬ 
cant.  The  most  important  issue  is  to  determine  the  effect  of  compressibility  on 
the  model  predicted  form  of  the  compressible  law  of  the  wall,  i.e.,  the  Karman 
constant,  k,  and  the  additive  constant,  B. 

In  addition  to  affecting  the  width  of  the  transition  region,  the  closure  coef¬ 
ficients  Rk  and  R^  determine  the  value  of  the  constant  in  the  law  of  the  wall, 
B,  for  a  turbulent  boundary  layer.  In  fact,  there  exist  unique  pairs  of  values 
{Rh,  Rw)  that  yield  B  =  5.0  for  the  incompressible  sublayer.  The  appropriate 
values  have  been  established  by  using  perturbation  methods  to  solve  the  sub¬ 
layer  equations  and  to  determine  B  from  the  limiting  form  of  the  solution  as 
y+  ->■  oo, 

Results  of  the  Mach  20  transition  application  above  suggest  that  larger  values 
of  Rjc  and  may  be  more  appropriate.  However,  there  is  a  possibility  that 
larger  values  may  be  appropriate  only  for  high  Mach  number  flows.  To  determine 
effects  of  compressibility,  we  have  redone  our  perturbation  analysis  of  the  viscous 
sublayer  with  effects  of  compressibility  included. 

For  coinpressible  flow,  the  sublayer  equations  assume  the  following  form. 

=  0  (6.1) 


d 

dy 


.  ^  du 


A 

dy  L 


,  .  du  (  u,  ut  \  dh  .  .  dk 


=  0 


d  \.  ^  .dk 

^\j.P  +  <rpT)- 


£ 

dy  L 


dy\ 

dy 


+  /^T 


.  V  dwl  f  du 

in  +  a-HTj  ~  +aa  p  I  — 


—  l3*puk  =  0 


•  jdpui  —  0 


where 


Pt  —  a*  pk/iii  and  ph  —  p^h^ 

We  can  integrate  Equations  (6.1  and  (6.2)  once,  wherefore 

,  .  du  2 

+  I^T)  =  PmU^ 


(6.2) 

(6.3) 

(6.4) 

(6.5) 

(6.6) 
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,  .  du  {  u  iiT  \  dh  ,  ^  dk 

ifi + j  ^  ^  (^-7) 

where  u,  is  friction  velocity,  is  density  at  the  surface,  and  is  surface 
heat  flux.  We  make  these  equations  dimensionless  by  introducing  the  following 
variables. 


+  «  +  UrV  ,+  Ic 

U  =  — >  r  =  - ,  k+  =  — 

Ur  Uu!  Uj 

+  p  h  ^  u^u 

P  =  1  n  =  — T,  Ijj  —  — 


f^W  /^t£> 


The  resulting  equations  are: 


^  =  1 


^[(;.+  +  ^>+)£^]+,4(^)  -r<.+«+i:+  =  o 


+ 


/?+*+  =  h+,  q+ 


(6.10) 


(6.11) 


(6.12) 


(6.13) 


(6.14) 


(6.15) 


This  is  a  sixth  order  system  of  coupled  nonlinear  equations,  and  we  thus 
need  six  boundary  conditions.  Four  conditions  can  be  applied  at  the  surface 
and  two  for  y"*"  — co.  Specifically,  from  no  slip,  the  velocity  and  turbulence 
kinetic  energy  vanish  at  the  surface  so  that 

,1  u+(0)  =  0  and  ;k+(0)  =  0  (6.16) 

1 

We  can  specify  either  surface  temperature  or  surface  heat  flux,  i.e., 

dh'^ 

/i+(0)  =  h+  or  ^^(0)  =  -^''LQw  (6.17) 

The  final  surface  boundary  condition  pertains  to  the  surface  value  of  the  specific 
dissipation  rate.  We  can  either  use  the  “rough-wall”  boundary  condition  for 
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which  the  surface  value  of  ui  is  finite  or  the  “smooth-wall”  condition  in  which 
the  asymptotic  limiting  form  for  w  is  imposed.  That  is  we  have: 


w+(0)=w+  or  ^  as  y+  ^0  (6.18) 

Boundary  conditions  far  from  the  surface  follow  from  matching  to  the  defect- 
layer  solution  or,  as  demonstrated  by  Wilcox, by  matching  to  the  wall  layer 
solution,  which  can  be  obtained  in  closed  form.  From  the  compressible  wall 
layer  solution  we  thus  conclude  that 


1  du'^ 


as  j/"*" 


oo 


(6.19) 


~  ^  2  ^  ^ ^  !/+-^oo  (6.20) 

where 

Mr  =  [{7-l)htr''"  (6.21) 

Finally,  the  velocity  in  the  wall  layer  is  given  by 

—  =  hny+  +  B  (6.22) 

Uj  K 

where  u*  is  the  scaled  velocity  that,  in  terms  of  the  dimensionless  variables  is 
given  by 


We  have  solved  Equations  (6.11)  through  (6.16)  subject  to  Equations  (6.16) 
through  (6.21)  by  modifying  our  incompressible  sublayer  program.  Computa¬ 
tions  have  been  done  with  the  primary  focus  on  the  effect  of  surface  temperature 
on  B,  which  is  obtained  from  the  solution  as  follows. 

'  r  1 

B  =  lim  - (.ny'^ 

^+—►00  K 

As  shown  in  Figure  6.1,  for  an  adiabatic  wall  and  with  R^,  Rk  and  R^,  given 
by  Equation  (2.16),  there  is  less  than  a  10%  change  in  B  for  varying  between 
500  and  10,000.  Although  more  computations  will  be  needed  to  consider  all 
possibilities,  it  appears  that  compressibility  has  a  relatively  small  effect  on  k—u 
model  predicted  sublayer  structure. 

Turning  to  the  ad  hoc  modification  discussed  in  Chapter  5,  the  constant  in 
the  law  of  the  wall,  B,  corresponding  to  using  Rk  and  R^  =  8.1  is  in  excess 


(6.24) 
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Figure  6.1:  Variation  of  the  constant  in  the  law  of  the  wall  with  snrface  tem¬ 
perature. 


of  10.  According  to  the  perturbation  solution,  the  value  of  B  can  be  reduced 
to  approximately  5.0  by  selecting  {Rk,Rw)  =  (18,28)  or  (/ZjtjJ?,,,)  =  (9.5, 8.1). 
Unfortunately,  both  choices  fail  to  yield  a  satisfactory  solution  for  the  Mach  20 
case  discussed  in  Chapter  5.  While  this  cursory  examination  of  the  problem  does 
not  exhaust  all  possible  choices  for  Rk  and  R^,  it  appears  that  some  modification 
of  the  damping  functions  for  a*  and  a  may  be  needed. 


6.2  The  Defect  Layer 

For  compressible  boundary  layers,  the  k  —  lj  model  equations  are  as  follows. 


d  ,  .  ^  .  n 


(6.25) 


du  du  dp  d  \,  .  dul 

dp  f  „  ,  d  \ f  u  PT  \  dk'\ 

^  -  Pt  I'p  +  (t*pt)^  (6.28) 

yy  \oyJ  ay  I  dy\ 


dh  dh  dp  ( _ 


(6.27) 


dk  dk  ( du 
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du  dw  W  a  2  ^  \r 

y.T  =  apk/u,  h  =  CpT,  P  -  pRT  (6.30) 

Viscous  eiFects  are  unimportant  in  the  defect  layer  so  that  the  closure  coefficients 
are  given  by 

=  1,  a  =  5/9,  (3*  =  9/100,  ^  =  3/40,  <7*  =  1/2,  <t  =  1/2  (6.31) 


To  analyze  the  compressible  defect  layer,  we  introduce  the  transformed  co¬ 
ordinates  ^  and  T)  defined  by 


-  ,,  I  r’  p  ,  .  v.K 

’'=amX  ft-'"’  ^w  =  — 


(6.32) 


where  Ut  and  pe  are  velocity  and  density  at  the  boundary- layer  edge,  Ut  is 
friction  velocity,  i  is  a  characteristic  length  scale  in  the  streamwise  direction, 
and  (S*  is  the  compressible  velocity  thickness  defined  by 


/■CO 

K=  - 

Jo  Pe 


(6.33) 


We  introduce  the  streamfunction,  ^{x,  y),  and  its  transformed  counterpart, 
T}),  defined  by: 

iP{x,y)  =  pMeAF{^,7,)  (6.34) 

Similaj-ly,  we  transform  the  various  flow  properties  according  to: 

u(x,y}  ~  UeU{^,Ti)-,  U^dF/dr]  ' 

PT{x,y)  =  p^UJ*N(^,tj) 


p{x,y)  =  PeRi^,1l) 
h{x,y)  = 


(6.35) 


The  transformed  equations  of  motion  possess  a  similarity  solution  in  which 
all  flow  properties  are  independent  of  x.  The  similarity  variable  is  y  and  the 
following  set  of  ordinary  differential  equations  define  the  similarity  solution. 

^  +  [ax  -  (2  -  M^)l3r  -  2^x1  =  I3tH  [l  -  U^]  (6.36) 
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WHK  (6.37) 

cr*4-  RN^{HK)  +[aT-{‘l-Ml)0T-^T]F4-{HK) 

dr]  Y  dr}  \  ^  di} 

=  [<j>T  -  -  4wt]  KHK  -  RN  -  WHK  (6.38) 

^^W^  dW 

fliV—  +[aT-{2-M^)PT-‘2uT\F~ 

=  [/?T-ar  +  4wT]^/l^- —  “f— "l  -4-^'*  (6-39) 

Pw  [  V«7-y  \a»?/  P 

and 

JV  =  a*-|  (6.40) 

where  is  Mach  number  at  the  boundary-layer  edge  and  the  four  parameters 
«T)  /?T)  wt  and  <i>T  are  defined  by  the  following  equations. 

2  dS:  ^  _KdP  _  K  dur  ^  _  26:  dT^ 

Cf  dx  ’  ^  Tu,  dx’  CfUr  dx  '  ^  CfTja  dx 

A  similarity  solution  exists  provided  these  four  parameters  are  independent 
of  X.  Although  we  have  not  yet  completed  our  analysis,  the  remaining  steps 
axe  straightforward  cuid  will  be  completed  during  the  next  year.  First,  we  must 
introduce  asymptotic  expansions  for  the  various  flow  properties.  Then,  exam¬ 
ination  of  the  momentum  integral  equations  should  tell  us  that,  similar  to  the 
incompressible  case,  Ur  varies  so  slowly  that  the  coefficient  uj-  is  negligible  to 
leading  order.  For  sifnilar  reasons,  <j)x  is  probably  negligible.  Further  manipula¬ 
tion  of  the  momentum  integral  equations  should  yield  a  relationship  between  ar 
and  /Jy.  The  resulting  equations  should  represent  the  similarity  solution  with 
the  equilibrium  parameter,  13^,  as  a  parameter.  The  equations  can  be  solved  nu¬ 
merically  with  straightforward  modifications  to  our  incompressible  defect-layer 
program. 


V' 
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Chapter  7 

Summary  and  Conclusions 


We  have  made  important  progress  in  developing  a  transition  model  during  the 
past  year.  Wilcox^®  presents  details  of  research  conducted  during  the  first  three 
months  of  the  project.  Briefly,  Wilcox^ describes  a  low-Reynolds-number  ver¬ 
sion  of  the  k  ~  ui  model  and  tests  it  for  fully- developed  turbulent  channel  and 
pipe  flow,  and  for  incompressible,  transitional  flat-plate  boundary-layer  flow. 

In  the  context  of  turbulent  boundary  layers,  this  report  shows  that  the  low- 
Reynolds-number  version  of  the  model  retains  all  of  the  best  features  of  the 
baseline,  high-Reynolds-number  version  of  the  model.  Most  notably,  the  model 
is  just  as  accurate  as  the  baseline  model  for  flows  with  adverse  pressure  gra^ 
dient.  This  guarantees  that  the  model  approaches  the  correct  post-transition, 
asymptotic  state. 

In  the  context  of  transitional  boundary  layers,  this  report  shows  that,  while 
offering  predictions  superior  to  those  of  simpler  correlations  designed  for  use 
with  algebraic  turbulence  models,  more  development  is  needed.  Most  impor¬ 
tantly,  transition  width  is  close  to  measured  width  only  for  low  transition 
Reynolds  numbers.  At  higher  transition  Reynolds  numbers,  model-predicted 
transition  is  much  more  abrupt  than  measured. 

While  an  ad  hoc  change  in  closure  coefficient  values  clearly  improves  matters 
for  a  Mach  20  case,  the  change  is  inconsistent  with  values  required  for  turbulent 
boundary  layers.  TMe  best  hope  for  resolution  of  the  problem  lies  in  the  per¬ 
turbation  analysis  of  Chapter  6.  Utilizing  the  sublayer  analysis,  we  should  be 
able  to  refine  the  viscous  damping  functions  to  arrive  at  suitable  low-Reynolds- 
number  limiting  forms  for  a*,  a,  /?*,  and  any  of  the  other  closure  coefficients. 

The  model’s  sensitivity  to  the  freestream  value  of  w,  to  some  extent,  miti¬ 
gates  the  model’s  insensitivity  to  initial  profiles.  The  defect-layer  analysis  offers 
promise  for  quantifying  and  eliminating  this  annoying  misfeature  of  the  model. 
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